Multimodal characterization of the visual network in Huntington's disease gene carriers

Objective A sensorimotor network structural phenotype predicted motor task performance in a previous study in Huntington’s disease (HD) gene carriers. We investigated in the visual network whether structure – function – behaviour relationship patterns, and the effects of the HD mutation, extended beyond the sensorimotor network. Methods We used multimodal visual network MRI structural measures (cortical thickness and white matter connectivity), plus visual evoked potentials and task performance (Map Search; Symbol Digit Modalities Test) in healthy controls and HD gene carriers. Results Using principal component (PC) analysis, we identified a structure – function relationship common to both groups. PC scores differed between groups indicating white matter disorganization (higher RD, lower FA) and slower, and more disperse, VEP signal transmission (higher VEP P100 latency and lower VEP P100 amplitude) in HD than controls while task performance was similar. Conclusions HD may be associated with reduced white matter organization and efficient visual network function but normal task performance. Significance These findings indicate that structure – function relationships in the visual network, and the effects of the HD mutation, share some commonalities with those in the sensorimotor network. However, implications for task performance differ between the two networks suggesting the influence of network specific factors

Emerging ensembles of kinetic parameters to identify experimentally observed phenotypes

Background: Determining the value of kinetic constants for a metabolic system in the exact physiological conditions is an extremely hard task. However, this kind of information is of pivotal relevance to effectively simulate a biological phenomenon as complex as metabolism. Results: To overcome this issue, we propose to investigate emerging properties of ensembles of sets of kinetic constants leading to the biological readout observed in different experimental conditions. To this aim, we exploit information retrievable from constraint-based analyses (i.e. metabolic flux distributions at steady state) with the goal to generate feasible values for kinetic constants exploiting the mass action law. The sets retrieved from the previous step will be used to parametrize a mechanistic model whose simulation will be performed to reconstruct the dynamics of the system (until reaching the metabolic steady state) for each experimental condition. Every parametrization that is in accordance with the expected metabolic phenotype is collected in an ensemble whose features are analyzed to determine the emergence of properties of a phenotype. In this work we apply the proposed approach to identify ensembles of kinetic parameters for five metabolic phenotypes of E. Coli, by analyzing five different experimental conditions associated with the ECC2comp model recently published by Hädicke and collaborators. Conclusions: Our results suggest that the parameter values of just few reactions are responsible for the emergence of a metabolic phenotype. Notably, in contrast with constraint-based approaches such as Flux Balance Analysis, the methodology used in this paper does not require to assume that metabolism is optimizing towards a specific goal

A critical evaluation of methods for the reconstruction of tissue-specific models

Under the framework of constraint based modeling, genome-scale metabolic models (GSMMs) have been used for several tasks, such as metabolic engineering and phenotype prediction. More recently, their application in health related research has spanned drug discovery, biomarker identification and host-pathogen interactions, targeting diseases such as cancer, Alzheimer, obesity or diabetes. In the last years, the development of novel techniques for genome sequencing and other high-throughput methods, together with advances in Bioinformatics, allowed the reconstruction of GSMMs for human cells. Considering the diversity of cell types and tissues present in the human body, it is imperative to develop tissue-specific metabolic models. Methods to automatically generate these models, based on generic human metabolic models and a plethora of omics data, have been proposed. However, their results have not yet been adequately and critically evaluated and compared. This work presents a survey of the most important tissue or cell type specific metabolic model reconstruction methods, which use literature, transcriptomics, proteomics and metabolomics data, together with a global template model. As a case study, we analyzed the consistency between several omics data sources and reconstructed distinct metabolic models of hepatocytes using different methods and data sources as inputs. The results show that omics data sources have a poor overlapping and, in some cases, are even contradictory. Additionally, the hepatocyte metabolic models generated are in many cases not able to perform metabolic functions known to be present in the liver tissue. We conclude that reliable methods for a priori omics data integration are required to support the reconstruction of complex models of human cells.Acknowledgments. S.C. thanks the FCT for the Ph.D. Grant SFRH/BD/ 80925/2011. The authors thank the FCT Strategic Project of UID/BIO/04469/2013 unit, the project RECI/BBB-EBI/0179/2012 (FCOMP-01-0124-FEDER-027462) and the project “BioInd - Biotechnology and Bioengineering for improved Industrial and Agro-Food processes”, REF. NORTE-07-0124-FEDER-000028 Co-funded by the Programa Operacional Regional do Norte (ON.2 - O Novo Norte), QREN, FEDER

Wanted dead or alive : high diversity of macroinvertebrates associated with living and ’dead’ Posidonia oceanica matte

The Mediterranean endemic seagrass Posidonia oceanica forms beds characterised by a dense leaf canopy and a thick root-rhizome ‘matte’. Death of P. oceanica shoots leads to exposure of the underlying matte, which can persist for many years, and is termed ‘dead’ matte. Traditionally, dead matte has been regarded as a degraded habitat. To test whether this assumption was true, the motile macroinvertebrates of adjacent living (with shoots) and dead (without shoots) matte of P. oceanica were sampled in four different plots located at the same depth (5–6 m) in Mellieha Bay, Malta (central Mediterranean). The total number of species and abundance were significantly higher (ANOVA; P<0.05 and P<0.01, respectively) in the dead matte than in living P. oceanica matte, despite the presence of the foliar canopy in the latter. Multivariate analysis (MDS) clearly showed two main groups of assemblages, corresponding to the two matte types. The amphipods Leptocheirus guttatus and Maera grossimana, and the polychaete Nereis rava contributed most to the dissimilarity between the two different matte types. Several unique properties of the dead matte contributing to the unexpected higher number of species and abundance of motile macroinvertebrates associated with this habitat are discussed. The findings have important implications for the conservation of bare P. oceanica matte, which has been generally viewed as a habitat of low ecological value.peer-reviewe

Multi-objective optimization of genome-scale metabolic models: the case of ethanol production

Ethanol is among the largest fermentation product used worldwide, accounting for more than 90% of all biofuel produced in the last decade. However current production methods of ethanol are unable to meet the requirements of increasing global demand, because of low yields on glucose sources. In this work, we present an in silico multi-objective optimization and analyses of eight genome-scale metabolic networks for the overproduction of ethanol within the engineered cell. We introduce MOME (multi-objective metabolic engineering) algorithm, that models both gene knockouts and enzymes up and down regulation using the Redirector framework. In a multi-step approach, MOME tackles the multi-objective optimization of biomass and ethanol production in the engineered strain; and performs genetic design and clustering analyses on the optimization results. We find in silico E. coli Pareto optimal strains with a knockout cost of 14 characterized by an ethanol production up to 19.74mmolgDW−1h−1 (+832.88% with respect to wild-type) and biomass production of 0.02h−1 (−98.06% ). The analyses on E. coli highlighted a single knockout strategy producing 16.49mmolgDW−1h−1 (+679.29% ) ethanol, with biomass equals to 0.23h−1 (−77.45% ). We also discuss results obtained by applying MOME to metabolic models of: (i) S. aureus; (ii) S. enterica; (iii) Y. pestis; (iv) S. cerevisiae; (v) C. reinhardtii; (vi) Y. lipolytica. We finally present a set of simulations in which constrains over essential genes and minimum allowable biomass were included. A bound over the maximum allowable biomass was also added, along with other settings representing rich media compositions. In the same conditions the maximum improvement in ethanol production is +195.24%

Enabling comparative modeling of closely related genomes: Example genus Brucella

For many scientific applications, it is highly desirable to be able to compare metabolic models of closely related genomes. In this short report, we attempt to raise awareness to the fact that taking annotated genomes from public repositories and using them for metabolic model reconstructions is far from being trivial due to annotation inconsistencies. We are proposing a protocol for comparative analysis of metabolic models on closely related genomes, using fifteen strains of genus Brucella, which contains pathogens of both humans and livestock. This study lead to the identification and subsequent correction of inconsistent annotations in the SEED database, as well as the identification of 31 biochemical reactions that are common to Brucella, which are not originally identified by automated metabolic reconstructions. We are currently implementing this protocol for improving automated annotations within the SEED database and these improvements have been propagated into PATRIC, Model-SEED, KBase and RAST. This method is an enabling step for the future creation of consistent annotation systems and high-quality model reconstructions that will support in predicting accurate phenotypes such as pathogenicity, media requirements or type of respiration.We thank Jean Jacques Letesson, Maite Iriarte, Stephan Kohler and David O'Callaghan for their input on improving specific annotations. This project has been funded by the United States National Institute of Allergy and Infectious Diseases, National Institutes of Health, Department of Health and Human Services, under Contract No. HHSN272200900040C, awarded to BW Sobral, and from the United States National Science Foundation under Grant MCB-1153357, awarded to CS Henry. J.P.F. acknowledges funding from [FRH/BD/70824/2010] of the FCT (Portuguese Foundation for Science and Technology) Ph.D. scholarship

Seagrass Canopy Photosynthetic Response Is a Function of Canopy Density and Light Environment: A Model for Amphibolis griffithii

A three-dimensional computer model of canopies of the seagrass Amphibolis griffithii was used to investigate the consequences of variations in canopy structure and benthic light environment on leaf-level photosynthetic saturation state. The model was constructed using empirical data of plant morphometrics from a previously conducted shading experiment and validated well to in-situ data on light attenuation in canopies of different densities. Using published values of the leaf-level saturating irradiance for photosynthesis, results show that the interaction of canopy density and canopy-scale photosynthetic response is complex and non-linear, due to the combination of self-shading and the non-linearity of photosynthesis versus irradiance (P-I) curves near saturating irradiance. Therefore studies of light limitation in seagrasses should consider variation in canopy structure and density. Based on empirical work, we propose a number of possible measures for canopy scale photosynthetic response that can be plotted to yield isoclines in the space of canopy density and light environment. These plots can be used to interpret the significance of canopy changes induced as a response to decreases in the benthic light environment: in some cases canopy thinning can lead to an equivalent leaf level light environment, in others physiological changes may also be required but these alone may be inadequate for canopy survival. By providing insight to these processes the methods developed here could be a valuable management tool for seagrass conservation during dredging or other coastal developments

MultiMetEval: comparative and multi-objective analysis of genome-scale metabolic models

Comparative metabolic modelling is emerging as a novel field, supported by the development of reliable and standardized approaches for constructing genome-scale metabolic models in high throughput. New software solutions are needed to allow efficient comparative analysis of multiple models in the context of multiple cellular objectives. Here, we present the user-friendly software framework Multi-Metabolic Evaluator (MultiMetEval), built upon SurreyFBA, which allows the user to compose collections of metabolic models that together can be subjected to flux balance analysis. Additionally, MultiMetEval implements functionalities for multi-objective analysis by calculating the Pareto front between two cellular objectives. Using a previously generated dataset of 38 actinobacterial genome-scale metabolic models, we show how these approaches can lead to exciting novel insights. Firstly, after incorporating several pathways for the biosynthesis of natural products into each of these models, comparative flux balance analysis predicted that species like Streptomyces that harbour the highest diversity of secondary metabolite biosynthetic gene clusters in their genomes do not necessarily have the metabolic network topology most suitable for compound overproduction. Secondly, multi-objective analysis of biomass production and natural product biosynthesis in these actinobacteria shows that the well-studied occurrence of discrete metabolic switches during the change of cellular objectives is inherent to their metabolic network architecture. Comparative and multi-objective modelling can lead to insights that could not be obtained by normal flux balance analyses. MultiMetEval provides a powerful platform that makes these analyses straightforward for biologists. Sources and binaries of MultiMetEval are freely available from https://github.com/PiotrZakrzewski/MetEv​al/downloads